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Abstract. High-mass microquasars are binary systems composed by a mas- 
sive star and a compact object from which relativistic jets are launched. Re- 
garding the companion star, observational evidence supports the idea that winds 
of hot stars are formed by clumps. Then, these inhomogeneities may interact 
with the jets producing a flaring activity. In the present contribution we study 
the interaction between a jet and a clump of the stellar wind in a high-mass 
microquasar. This interaction produces a shock in the jet, where particles may 
be accelerated up to relativistic energies. We calculate the spectral energy distri- 
butions of the dominant non-thermal processes: synchrotron radiation, inverse 
Compton scattering, and proton-proton collisions. Significant levels of X- and 
7-ray emission are predicted, with luminosities in the different domains up to 
~ 10 34 - 10 35 erg s 1 on a timescale of about ~ 1 h. Finally, jet-clump inter- 
actions in high-mass microquasars could be detectable at high energies. These 
phenomena may be behind the fast TeV variability found in some high-mass X- 
ray binary systems, such as Cygnus X-l, LS 5039 and LS 1+61 303. In addition, 
our model can help to derive information on the properties of jets and clumpy 
winds. 



1. Introduction 

The mass loss in massive stars is thought to take place via supersonic inhomo- 
geneous winds. Considerable observational evidence supports the idea that the 
wind structure is clumpy (e.g. Owocki & Cohen 2006), although the proper- 
ties of clumps are not well-known as a cosequense of the high spatial resolution 
necessary for the clumps detection. Some massive stars are accompanied by a 
compact object and present transfer of matter to it forming an accretion disk 
and, in high-mass microquasars (HMMQs), bipolar relativistic outflows. 

Non-thermal emission has been observed in microquasar jets from radio to 
X-rays, and it is thought that 7-rays could be also produced in jets. Recently, a 
TeV flare has been detected from the HMMQ Cygnus X-l (Albert et al. 2007). 
Transient 7-ray events could have been detected as well from the high-mass 
X-ray binaries LS 5039 (Aharonian et al. 2005), LS 1+61 303 (Albert et al. 
2006) and perhaps from the HMMQ Cygnus X-3 (ATels 1492, 1547 and 1585). 
Some authors have been suggested that this strongly variable 7-ray emission 
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can be produced by the interaction between the jet and the stellar wind of the 
companion star (e.g. Owocki et al. 2009, Araudo et al. 2009). 

In the present contribution we propose a model to explain these 7-ray 
flares based on the interaction between the jets of a HMMQ with wind inhomo- 
geneities. The clumps can eventually penetrate in the jet, leading to transient 
non-thermal activity that may release a significant fraction of the jet kinetic lu- 
minosity in the form of synchrotron, inverse Compton (IC), and proton-proton 
(pp) emission. 

2. The scenario 

To study the interaction between a clump of the stellar wind and a jet in a 
HMMQ, we adopt a scenario with similar characteristics to Cygnus X-l. We 
fix the separation between the compact object and the massive star to a = 
3 x 10 12 cm. The stellar mass loss rate is adopted to be M+ = 3 x 10 -6 M@ yr _1 , 
with a terminal wind velocity v w ~ 2.5 x 10 8 cm s . A sketch of the scenario 
is presented in Figure [TJ 

2.1. Clump and jet model 

The clump is taken to be spherical and we consider two values of it radius: 
R c = 10 10 and 10 11 cm. We have assumed that the velocity of the clumps is 
equal to the velocity of the wind, i.e v c = v w . In order to obtain a large density 
contrast between the clump and the jet we assume dense and homogeneous 
clumps with n c = 10 12 cm -3 , which corresponds to a clumpy wind filling factor 
of / = M*/47ra 2 mp'U c n c = 0.005, where m p is the mass of the proton. 

We consider a jet dynamically dominated by cold protons with a bulk ve- 
locity Vj = 0.3 c and a kinetic luminosity L^in ~ 3 x 10 36 erg s _1 (e.g. Gallo et 
al. 2005). We assume here that the jet radius is one tenth of the jet height, i.e. 
Rj(z) = 0.1 z, and the jet-clump interaction is taken to occur at z- m t = a/2. The 
density of the jet material at z- mt results nj = 4.7 x 10 7 cm -3 and then, the ratio 
between the clump and the jet densities is x = 2.1 x 10 4 . This parameter will 
be very relevant for the jet-clump interaction estimates. 

2.2. Dynamics of the interaction 

In order to study the physical processes of the jet-clump interaction, we consider 
the collision of a single clump with the jet. The huge density contrast x allows 
the clump to cross the boundary of the jet and fully penetrate into it. In the 
context of this work, thermal conduction, clump expansion, magnetic fields and 
gravitational forces are not dynamically relevant for the jet-clump interaction 
and will be neglected. 

In the Laboratory Reference Frame (LRF), the clump will take a time t c to 
fully enter the jet and a time ij to cross it roughly at the wind velocity: 

t c ~ 2R c /v c ~ 80 s (Rc = 10 10 cm) and ~ 800 s (R c = 10 11 cm), (1) 

2 R' 

U 1 = 1.2 x 10 3 s. (2) 

Vc 
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When the clump interacts with the jet, a shock is formed in the clump 
and propagates in the direction of the jet motion with a velocity v cs ~ v \l \/X- 
The clump-crossing time, which is the characteristic timescale of the jet-clump 
interaction, can be defined as: 

2 R 

t cc ~ 3 x 10 2 s (R c = 10 10 cm) and ~ 3 x 10 3 s (R c = 10 11 cm). (3) 

v cs 

A shock (the bow shock) is also formed in the jet when its material collides 
with the clump. We assume the clump/bow-shock separation distance as x ~ 
0.2R C (van Dyke & Gordon 1959), and thus the time required to reach the steady 
state regime results: 

0.2^ 0.2R C 5 t cc 
tbs jj- ~ « t cc . 4 



Once the clump is inside the jet, the latter accelerates the clump up to the 



background velocity, Vj, with an acceleration g ~ v?/xRc (e.g. Fragile et al. 
2004) in a time 

V' 

tg ~ — ~ y/X t cc » t cc - (5) 

When one fluid exerts a force against another fluid of different density, the 
hydrodynamical Rayleigh- Taylor (RT) instability eventually develops, leading 
to the perturbation and potential disruption of the clump. On the other hand, 
after the bow shock is formed, the jet material surrounds the clump and we have 
two fluids with a large relative velocity (in our case ~ v\). This situation leads 
to Kelvin-Helmholtz (KH) instabilities. The timescale for these instabilities, 
considering that the length of the perturbed region is ~ R c , are: 



tRT ~ \jRc/g ~ t cc and t KH ~ R ° — ~ t cc . (6) 

v j 

We assume that for a t c ~ t cc /5 we can treat the jet-clump interaction 
as the clump were fully inside the jet. At this stage we do not consider the 
penetration of the clump into the jet, and assume in our treatment that the 
former is completely inside the latter (i.e. the system has cylindrical symmetry; 
see Fig. [1]). We notice that the clump might not escape the jet if tj > £rt/kh 
due to disruption through instabilities, although numerical simulations show 
that the instability timescales are actually longer than the clump crossing time 
by a factor of a few ( e.g. Klein, McKee & Colella 1994). 



3. Non-thermal processes 

Given the characteristics of our scenario, the shock in the clump is radiative, 
and the bow shock is adiabatic. For these reasons, the shocked and heated 
material of the clump will radiate (by thermal Bremsstrahlung) a non-negligible 
part of the energy tranferred by the jet, whereas in the bow shock non-thermal 
processes take place. We focus here in the non-thermal emission. 
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Compact Obji 



Figure 1. Sketch of a HMMQ with a clumpy wind (left) and the jet-clump 
interaction (right). 



3.1. Particle acceleration and radiative cooling 

Assuming that non-relativistic diffusive (Fermi I) shock acceleration takes place 
in the bow-shock, electrons and protons with a charge q will be accelerated up 
to an energy E e ^ p in a time 




'—l-l lite* (?) 

in the Bohm limit and for perpendicular shocks (Protheroe 1999). The magnetic 
field in the acceleration region is B\^ s , and we will consider two values for it. First, 
we consider a value obtained assuming that the magnetic energy density is a 10 % 
of the plasma internal energy density downstream and we obtain Bb s = 150 G. 
We adopt a sub-equipartition value to make B^ dynamically negligible with 
respect to matter. In addition, we adopt B^ s = 1 G to check the impact on our 
results of a magnetic field much smaller than in the sub-equipartition case. 

In the scenario presented in this work, relativistic leptons lose their en- 
ergy mainly by synchrotron radiation and IC scattering. The former have the 
following cooling time 

4.1 x 10 2 

At Z; n t, the energy density of the photons from the star is itph — 2.4 x 10 erg 
cm -3 , being the typical photon energy eo ~ 10 eV. For y = eoE e /(5.1 x 
10 5 eV) 2 > 1 the IC interaction takes place in the Klein-Nishina (KN) regime. 
A formula for the IC cooling time valid in both Thompson (Th) and KN regime 
under a photon field with a narrow energy distribution is (e.g. Bosch-Ramon & 
Khangulyan 2009): 

= 6.1 x 10 12 6q (1 + 8.32/) (1 + lV) 
IC u ph ln(l + 0.2y) (l + 0.5y + 1.3y 2 ) S ' U 
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Regarding hadronic emission, 7-rays are produced if relativistic protons 
interact with nuclei through inelastic pp collisions, beeing the proton cooling 
time 

2.2 x 10 15 . , 

hp = ~ J—_ / — K _ \\ s - ( 10 ) 



4n; (o^ + Omini-r^- 

J \ \l.lm p c 



In the bow-shock region, this process is negligible. Otherwise, if relativistic 
protons accelerated in the bow shock penetrate into the clump, pp interactions 
can become an efficient process to generate 7-rays because n c 3> nj. 

Besides radiative cooling, the escape of particles from the acceleration region 
is another kind of energy losses. The corresponding timescale, r esc , takes into 
account the advection of relativistic particles by downstream bow-shock material 
(^adv ~ -Rc/^jps) and the diffusion of particles (tdiflf ~ x 2 /2Db, where Db is the 
Bohm diffusion coefficient), being r esc = minjtdifr, t a d v }- Equating energy gains 
and losses, the maximum energy achieved by electrons (-E™ ax ) accelerated in the 
bow shock is estimated (see Fig. [2]). On the other hand, the maximum energy 
for protons (E™ x ) is constrained by the size of the acceleration region. In Table 
[U £™ ax and E™ ax are shown. 




10 12 

log(E, [eV]) 




Figure 2. Acceleration and radiative loss time for electrons in the bow- 
shock region. The advection and diffusion times are shown for R c = 10 10 
(thin line) and 10 11 cm (thick line). 



Table 1. Maximum energies achieved by acelerated particles in the bow 
shock and the energy break of electrons. 



R c [cm] 


10 iU 


10 iU 


10 n 


10 n 


B hs [G] 


1 


150 


1 


150 


^max [ eV ] 


1.5 x 10 11 


8 x 10 11 


1.5 x 10 iZ 


8 x 10 11 


Eg** [eV] 


6 x 10 11 


9 x 10 13 


6 x 10 12 


9 x 10 14 


E h [eV] 


4.5 x 10 10 


3 x 10 9 


4.5 x 10 11 


3 x 10 8 
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3.2. High energy emission 

In this section we calculate the Spectral Energy Distribution (SED) of the emis- 
sion produced by the most relevant non-thermal radiative processes: synchrotron 
and IC radiation in the bow-shock region and pp in the clump. 

We assume an injected population of relativistic particles in the bow-shock 
region that follows a power-law energy distribution of the form: 

Qe, P ( E e,p) = K e , P E~ 2 p exp(-£ e , p / E^) , (11) 

typical for Fermi I acceleration. The normalization constant K e>p is determined 
assuming that the (non-thermal) luminosity of accelerated particles is Lnt = 
0.25(i? c /i?j) 2 Li c i n (note that the predicted fluxes scale linearly with the adopted 
non-thermal fraction 0.25). To estimate the particle energy distribution N ep , 
we solve the kinetic equation in the one-zone model approximation for the bow- 
shock region (e.g. Ginzburg & Syrovatskii 1964): 

Ot dE £tP dt * T esc 

The relativistic leptons reach the steady state well before the shock has 
crossed the clump. Particles downstream with low energies escape advected in 
the shocked material of the jet before cooling radiatively, producing a break in 
the energy spectrum of particles. From i a d v = t syn and t a d v = ^diff we determine 
the electron break energies E^. The values are shown in Table [TJ The most 
energetic electrons can diffuse up to the clump (B^ a = 1 G) or lose their energy 
inside the bow-shock region by synchrotron and IC radiation (B^ s = 150 G). We 
note that the electrons with energies E e > E^ can reach the clump and radiate 
inside it. In this work, we do not estimate the lepton emission in the clump. 

As noted in Sect. 13.1.1 relativistic protons do not suffer significant pp losses 
in the bow-shock region (i^is <C t pp ). These protons can also reach the clump if 
they are not advected by the shocked material of the jet. Imposing that t^iflf < 
i a d v , the minimum energy necessary to reach the clump is -E™ in = 0.025£™ ax , 
and the maximum energies of these protons are shown in Table [TJ On the 
other hand, to confine relativistic protons in the clump requires a very large 
magnetic field and we assume that protons will cross the whole clump in a time 
~ R c /c < tpp. Then, the distribution of relativistic protons in the clump is 
N p {Ep) = (R c /c)Q p (E p ). 

With the steady distributions of relativistic electrons in the bow shock, 
N e (E e ), and protons in the clump, N p (E p ), we calculate the SEDs. We estimate 
the specific luminosity as 

eL(e) = e / N e , p {E e , p ) j(E e , p , e) dE e , p , (13) 

e,p 

where j(E e ,e) is the emissivity of the process and e is the photon energy. 

Synchrotron and IC emissivities are computed using the standard formulae 
given in Blumenthal & Gould (1970). For the later process we consider the 
parametrized differential cross section valid in both Th and KN regimes and 7r°- 
decay emissivity is calculated following Kelner et al. (2006). The synchrotron 



High-energy Bares from jet-clump interactions 



7 



emission is computed in the optically thin case, neglecting the impact of syn- 
chrotron self-absorption, but at high energies (IC and pp) 7 — 7 absorption is 
taken into account. 

Synchrotron emission is dominant at X-rays for B^ = 150 G, with L syn ~ 
10 35 erg s _1 . Inverse Compton scattering produces 7-rays up to very high ener- 
gies, dominating the radiation output in the cases with B^ = 1 G. In our calcu- 
lations, the highest luminosity achieved is L\q ~ 10 35 erg s _1 for R c = 10 11 cm, 
although 7 — 7 absorption can reduce substantially the emission above 100 GeV. 
Proton-proton collisions in the clump can also produce 7-rays at energies that 
may be as high as ~ 10 14 eV (B^ s = 150 G). The maximum luminosity ob- 
tained by pp is nevertheless quite modest, L pp ~ 10 32 erg s _1 for R c = 10 11 cm, 
although denser and/or bigger clumps, and more powerful jets may yield de- 
tectable amounts of photons outside the 7 — 7 absorption range (0.1-10 TeV). 
The computed SEDs are shown in Fig. [3J only for the cases with B^ s = 150 G. 



T 




J 1 I 1 I 1 I 1 I 1 I 1 — 1 1 I 1 I 1 I I I I , I , I , I 1 I , I , I , I , I , L 

-4 -2 2 4 6 8 101214 -4 -2 2 4 6 8 101214 

l°9(e [eV]) |og(e [eV]) 



Figure 3. Leptonic and hadronic emission from the bow shock and from 
the clump, respectively, for Bb s = 150 G, and R c = 10 10 (left) and 10 11 cm 
(right). The curves of both absorbed (thick lines) and unabsorbed (thin lines) 
IC and pp radiation are shown. 



4. Effect of many clumps 

In general, many clumps can interact with the jet and is possible to have several 
clumps inside it simultaneously (see Owocki et al. 2009). However, at z < z c = 
4 x 10 11 cm the destruction time t cc is less than the penetration one (t c ) and 
the clump will be destroyed before completely entering into the jet. Then, the 
model presented here is not valid at z < z c (neither if Rj < R c ). 

Considering a canonical geometry for the jet and a filling factor of the wind 
/ = 0.005, we estimate the number of clumps -N c - inside the jet, from z = z c 
up to z ~ a. We obtain iV c ~ 350 and ~ 0.5 for R c = 10 10 and 10 11 cm, 
respectively. As a consecuence, flares produced by jet-clump interactions could 
be a sporadic phenomenon, for a low N c (R c = 10 11 cm), or may appear as a 
modulated steady activity, for a high iV c (R c = 10 10 cm). In the latter case, 
the resulting SED will be similar than the obtained in the previous section, but 
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multiplied by iV c . Nevertheless, we note that the jet may be disrupted in those 
cases when too many clumps are simultaneously present inside the jet (Araudo et 
al. 2009). However, more detailed calculations of the dynamics of the jet-clump 
interaction are required to clarify this issue. 

5. Conclusions 

In the present work, we study the main physical processes and the radiative 
outcomes from the interaction between the jet of a HMMQ with a clump of the 
companion stellar wind. This interaction produces a bow-shock in the jet where 
particles are accelerated up to very- high energies. The relativistic electrons are 
cooled efficiently by synchrotron and IC radiation in the bow shock, and protons 
by 7r°-decay from pp interactions in the clump. 

Due to characteristics of the interaction, the expected emission can be tran- 
sient if N c is low. Given the typical dynamical timescales, the flare duration may 
last from minutes to a few hours. The flares of jet-clump interactions would have 
associated lower (synchrotron) and higher energy components (IC, pp emission). 
The total level of emission, the importance of the different components, and the 
duration of the flares, can give information on the jet power and size, clump size 
and density, and magnetic fields in the interaction regions. Therefore, besides 
the jet itself, the clump properties can be probed by observations at high en- 
ergies of transient activity in HMMQs, being a new tool to study the winds of 
massive stars. 
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